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OS ■ Abstract 
0^ ■ 

, Numerical data |l^ for the SU{2) confining string in the maximal abelian pro- 

I jection are analysed. The distribution of the electric flux and monopole currents 

On ■ are perfectly described by the classical equations of motion for the dual Abelian 

Higgs model. The mass of the vector boson is equal to the mass of the monopole 
I (Higgs particle) within numerical errors. The classical energy per unit length of the 

Abrikosov vortex reproduces 94% of the full non-Abelian string tension. 



1. The monopole confinement mechanism [|I| in lattice gluodynamics is confirmed by 
many lattice numerical calculations (see reviews and references therein). The simplest 
5^ I effective dynamical picture corresponding to this mechanism is the dual Abelian Higgs 
model (AHM) with the Lagrangian: 

^AHM = -^A{B) + \\{d, - ^B,W + m? - ri'f , (1) 

where is the dual gauge field interacting with the monopole field Electrically charged 
"quarks" can be introduced using the Zwanziger formalism Q]. The confinement of 
electric charges is due to the dual Abrikosov string and exists already at the classical level. 
The first fact in favor of an effective Lagrangian like (|l]) was the numerical proof of the 
existence of the monopole condensate in the confining phase of lattice gluodynamics . 
The further numerical investigation allowed to derive an effective monopole Lagrangian 
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which is in agreement with (|I|). There are various attempts to determine the parameters 
of the effective dual theory from phenomenological analysis and from calculations on 
the lattice (e.g. 0, ^). The results can be expressed through the gauge boson mass 
my = gf] and Higgs (monopole) mass tuh = 2y/2Xri'^. In most cases, the analyses have 
led to the result my ~ ttlh- This is the so-called Bogomolnyi limit of the AHM, where 
several exact results can be obtained. In particular the classical Abrikosov string tension 
is: 



a=J-^m?, m = rriH = my . (2) 

V r 

Of course the simple model (|ID cannot describe all properties of QCD, and it is im- 
portant to find the region of applicability of this effective theory. An interesting example 
was found by Bali, Schlichter and Schilling |]10| who have presented results on the profile 
of the flux tube between heavy quarks in SU{2) lattice gauge theory. Working within 
the maximal Abelian projection on a 32^ lattice at /5 = 2.5115, they have studied the 
correlators of the Abelian electric field E and of the density of Abelian monopole cur- 
rents k with Abelian Wilson loops. Such correlators, divided by the expectation values 
of the Wilson loops themselves, describe the electric field and monopole currents induced 
by the presence of a static chromoelectric string. For Wilson loop in the zt plane only 
the z component Ez{p) of E and the azimuthal component ke{p) of k were found to be 
non-vanishing (here p = y/x'^ + y'^ is the distance from the center of the flux tube and the 



azimuthal angle 6 is defined as usual tan 6' = y/x). The authors of [T^ fit the numeri- 
cal data using specific parameterization of Ez{p) and kg{p) which respect the boundary 
conditions for the continuum dual Abrikosov string. 

2. Below we continue the analytical analysis of the numerical data for the profile of 
the confining string. We will avoid to use specific ansatze. Instead, we work directly with 
the Ginzburg-Landau equations, fitting the parameters of the dual AHM in such a way 
that the dual Abrikosov string reproduces the numerical data for Ez{p). Our experiments 
revealed that the usual l^l'' form of the Higgs field potential (|l]) is sufficient. The inclusion 
of higher order term Cgl^l^ does not change the results and gives cg = 0. In the cylindrical 
coordinate system t, z, p, p^ = XaXa {a = 1, 2), for the unitary gauge 53m $ = 0, the static 
axially symmetric and infinitely long ANO vortex can be parameterized by two functions 
b{p) and /(p): 



Ba = Sab J ^- - bip)j , $ = vfip) , Bo = B, = 0. (3) 
The energy per unit length [i.e. the string tension a) can be easily derived from (^: 

cr = ^Jpdp\ li-ipb]'? + iKUr + \Kn- - hf + \mWAf - I? (4) 

9 \ P P ) 
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The variational principle gives rise to the following equations of motion 



l.[pb']'-^b + mlf(-^-b]=0, (5) 



^[p/'r = /(^-&)' + ^</(f -1), (6) 
which should be supplied with the boundary conditions 

6(0) = /(0)=0, lim b{p) = -^, jim/(p) = l. (7) 

Note that (j^) is nothing but the dual Ampere law, k = curl^, in cylindrical coor- 
dinates. We have numerically solved the classical equations (|^J^) fitting the data for 
Ez{p) of Ref. II^. The fitting curves are plotted in Fig. 1. The numerical values of the 
parameters of the fit are the following : 

g/2Ti = 0.9274 ±0.0066, 

mv = 0.5733 ±0.0383 = (1.3123 ±0.0771) GeV, (8) 
mu = 0.5945 ± 0.0063 = (1.3614 ± 0.0143) GeV . 

The vector and scalar masses are equal to each other within the numerical errors. The 
physical meaning of this equality deserves a separate discussion although, as we mentioned, 
several models predict it. From a practical point of view we can estimate the string tension 
using eq. (|^): 

= 0.1744 ±0.0231 ^ (400.1 ±53.0) MeV ^ 0.91 ^/asu (2) ■ (9) 

Here and in (H) the dimensionless quantities correspond to the lattice spacing a = 1, to 
get the dimensional quantities we use the lattice results |jTO[ corresponding to j3 = 2.5115: 
a = 0.086 fm, ^su{2) = 440 MeV. 

3. One sees from Fig. 1 that, in spite of the good approximation to Ez{p), the 
magnetic current distribution is not perfectly well reproduced by our model in the region 
1 < p < 3. Since the Ampere law (^ is an integral part of the dual AHM, we conclude that 
either the dual description is invalid in this region or the numerical data are not reliable 
here due to finite lattice spacing effects. Ideally, to check the validity of dual AHM 
one should perform the same lattice measurements for larger values of (3 when lattice 
configurations are closer to the continuum. Such an enterprise is certainly worthwhile 
although is impossible practically at the moment. Instead, we try to estimate the effect 
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of hypercubic discretization by repeating our procedure on a coarse lattice similar to the 
one used in Ref. The discretized action corresponding to (|l]) is 

^ = 7^ E Fl,. + ^ E 1^^ - ^^-f^r + A Ed^^l' - v'f , (10) 

y x,fj.u x,fi X 

where F^^j^^i, = B^.^^ + B^^f^^y — B^+o^^ — B^^,, and a = 1. Since we consider an infinitely 
long ANO vortex the problem becomes essentially two-dimensional. In the unitary gauge 
^e^x = vfx, '^fn^x = 0, the string tension is 

^ = A(e(W + ^v E lU-e'""-^ fx^,\' + -mWvT.if!-^A^ (11) 

"^9 y X x,At=l,2 ^ X J 

Contrary to the continuum considerations above there are no apriori ansatze for the 
ANO string solution on the lattice. Nevertheless, it is known how to introduce a single 
vortex to the system |TT|: the recipe is to make a shift -Fx,i2 -^z,i2 + '^t^Sx,xo 
action (|TT|) . Of course, periodic boundary conditions are not allowed anymore. We have 
solved the corresponding equations of motion on the square lattices 32^ -=-70^ and with free 
boundary conditions. We did not find any finite volume effects. Note that a 32^ lattice 
in our analysis corresponds to the 32^ lattice of ref.|l^. We found that the numerical 



data for Ez and kg are strikingly well reproduced by a classical solution of the model 
The discrepancy with the observed magnetic currents distribution mentioned above 
disappears completely. Fig. 2 represents the results of the fit. The corresponding values 
of the parameters are: 

g/27r = 0.9519 ±0.0041, 

mv = 0.4522 ±0.0206 = (1.0351 ±0.0472) GeV, (12) 
tuh = 0.4747 ± 0.0600 = (1.0866 ± 0.1373) GeV . 

It appears that vector and scalar masses are equal to each other again although they are 
now ~23% smaller than in (^. Thus the seeming discrepancy between SU{2) gluody- 
namics and classical effective dual AHM is recognized as merely a discretization artifact. 
The string tension (0) turns out to be 

^/a = 0.1808 ±0.0213 ^ (414.8 ± 48.9) MeV ^ 0M^/asui2) ■ (13) 

4. To summarize, we have shown that the Abelian electric field profile and the mag- 
netic current distribution around the static confining string in SU (2) gluodynamics may 
be well described as a classical ANO vortex. Moreover, we found that the classical energy 
per unit length of the ANO vortex accounts for 94 % of the quantum SU (2) string tension. 
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The effective dual Abelian Higgs Model appears to lie on the border between type-I and 
type-II superconductivity. The magnetic charge of the Higgs field (|T^) corresponds to an 
electric charge of the quark q 0.95 (the charge of the quark in SU{2) Wilson loop is 
unity). This approximate correspondence of classical AHM and quantum SU{2) gluody- 
namics is dynamically more specific than just Abelian dominance [Q. The detailed study 



of this correspondence implies investigation of the Abrikosov vortex in quantum AHM, 
which is in progress now. 
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Figure Captions. 

Fig. 1 The fit of the data of Ref. |T0[ by a continuum solution. 

Fig. 2 The same as in Fig. 1 but using a classical solution on the two-dimensional lattice. 
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